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ABSTRACT 

We  present  global  potential  energy  surfaces  for  the  three  lowest  triplet  states  in  0(  P)  + 
H20(X  Ai)  collisions  and  present  results  of  classical  dynamics  calculations  on  the  0(  P)  + 
H20(X  A  i  )  — >  OH(X  II)  +  OH(XTl)  reaction  using  these  surfaces.  The  surfaces  are  spline-based 
fits  of  -20,000  fixed  geometry  ab-initio  calculations  at  the  CASSCF+MP2  level  with  a 
0(4s3p2dlf)/H(3s2p)  one  electron  basis  set.  Computed  rate  constants  compare  well  to 
measurements  in  the  1,000-2,500  K  range  using  these  surfaces.  We  also  compute  the  total,  ro- 
vibrationally  resolved,  and  differential  angular  cross  sections  at  fixed  collision  velocities  from 
near  threshold  at  -4  km  s'1  (16.9  kcal  mol"1  collision  energy)  to  11  km  s"1  (122.5  kcal  mol"1 
collision  energy),  and  we  compare  these  computed  cross  sections  to  available  space-based  and 
laboratory  data.  A  major  finding  of  the  present  work  is  that  above  -40  kcal  mol"1  collision 
energy  ro-vibrationally  excited  OH(X“IT)  products  are  a  significant  and  perhaps  dominant 
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contributor  to  the  observed  1-5  p  spectral  emission  from  0(  P)  +  tBCXX  Ai)  collisions.  Another 
important  result  is  that  OH(X  II)  products  are  fonned  in  two  distinct  ro-vibrational  distributions. 
The  ‘active’  OH  products  are  formed  with  the  reagent  O-atom,  and  their  ro-vibrational 
distributions  are  extremely  hot.  The  remaining  ‘spectator’  OH  is  relatively  ro-vibrationally  cold. 
For  the  active  OH,  rotational  energy  is  dominant  at  all  collision  velocities,  but  the  opposite  holds 
for  the  spectator  OH.  Summed  over  both  OH  products,  below  -50  kcal  mol1  collision  energy, 
vibration  dominates  the  OH  internal  energy,  and  above  -50  kcal  mol"1  rotation  is  greater  than 
vibrational  energy.  As  the  collision  energy  increases,  energy  is  diverted  from  vibration  to  mostly 
translational  energy.  We  note  that  the  present  fitted  surfaces  can  also  be  used  to  investigate  direct 
collisional  excitation  of  H20(X*Ai)  by  0(3P)  and  also  OH(X2n)  +  OH(X2n)  collisions. 
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I.  INTRODUCTION 
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Collisions  of  0(  P)  with  H20(X,  Ai)  give  rise  to  spectral  radiation  observed  from  the  long 
wave  infrared  (-20-10  p)1-3,  through  the  infrared  (10-1  p)1-9  and  up  to  the  vacuum  ultraviolet 
(-0.25  p), 10-12  depending  on  the  relative  collision  energy.  These  emissions  are  important  to 
characterize  because  they  can  be  a  significant  background  noise  source  for  telescopes  mounted 
on  spacecraft  in  low-Earth-orbit  or  LEO.  In  LEO  (-300  km  altitude  and  -8  km  s-1  relative 
atmospheric  velocity),  0(3P),  which  is  the  main  component  of  the  atmosphere  at  these  altitudes, 
collides  with  H20(X,1Ai)  from  outgassing  spacecraft  surfaces  or  from  spacecraft  engine 
exhaust.4 
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The  mechanisms  that  give  rise  to  the  observed  radiation  in  0(  P)  +  H20(X,  Ai)  collisions 
remain  elusive,  however,  because  the  rarefied  conditions  and  large  relative  collision  velocities 
necessary  are  difficult  to  achieve  in  the  laboratory,  and  so  experimental  data  remain  sparse  and 
difficult  to  interpret.  Focusing  on  the  1-5  p  region,  Bernstein  et  al.7  have  shown  through  analysis 
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of  space  experiments  and  laboratory  data  that  the  observed  radiation  from  0(  P)  +  H20(X,1A1) 
collisions  may  be  consistent  with  forming  ro-vibrationally  excited  OH  through  hydrogen 
abstraction, 

0(3P)  +  H20(X,1A1)  ->  OH(X,2n)  +  OH(X,2n)  ,  AH  =  +16.9  kcal  mol'1,  (1) 

and/or  fonning  ro-vibrationally  excited  H2O  through  direct  collisional  excitation, 

0(3P)  +  H20(X,1A1,  (v1v2v3,JK))  ->  0(3P)  +  H20(X,1A1(v1  v2  v3’,JK )),  (2) 

depending  on  the  collision  velocity.  Other  studies6  have  attributed  all  such  radiation  from  1-5 
p  to  Reaction  (2).  It  is  particularly  difficult  to  distinguish  between  (1)  and  (2)  in  the  1-5 
p  spectral  region  because  spectral  observations  remain  the  main  source  of  data  and  0H(X,21T) 
and  FEOtX,1  A|)  have  highly  overlapping  spectral  distributions  in  this  region.  Furthermore,  since 
spectral  observations  are  only  sensitive  to  internally  excited  products,  the  magnitude  and  energy 
dependence  of  these  observations  minimize  contributions  from  products  formed  with  little 
internal  excitation,  products  which  may  in  fact  carry  most  of  the  total  reaction  cross  section  and 
dynamical  information.  In  addition,  theoretical  examination  of  the  reactions  (1)  and  (2)  has  been 

2 

"Approved  for  public  release;  distribution  unlimited." 


problematic  because  they  involve  computation  of  high-spin  open  shell  wave  functions  which 
presents  special  difficulties  for  theory.  There  have  been  theoretical  studies  which  used  non- 
reactive  0(  P)  +  fbCXX,  Ai)  potential  surfaces  to  examine  the  collisional  excitation  channel  (Eq. 
2).  '  Although  these  studies  have  had  some  success  in  explaining  measurements  in  the  1-5  p 
region,  there  remain  major  differences  between  these  studies  and  measurements  at  certain 
collision  velocities  and  significant  differences  between  the  theoretical  results  themselves. 

Besides  being  important  in  LEO  conditions,  collisions  of  0(  P)  +  H20(X,  Ai)  are  also 
important  in  H2/O2  combustion  chemistry.17  The  reverse  of  Reaction  (1)  is  an  important  step  in 
such  chemistry.  In  fact,  the  surfaces  involved  in  Reaction  (1)  are  part  of  a  family  of  triplet  H2O2 
surfaces  leading  to  H2(XIIA)  +  02(X3Zg')  and  H(2S)  +  OOH(X2A’),  pathways  critical  in 
understanding  such  combustion  chemistry.  The  potential  surfaces  involved  in  0(  P)  + 
H20(X,1Ai)  collisions  have  been  probed  recently  using  transition  state  photoelectron 
spectroscopy,18’19  which  has  given  information  about  the  barrier  height  and  other  transition  state 
properties  for  Reaction  (1).  Instead  of  translational  energy  as  in  previous  space  and  laboratory 
measurements,  Ref.  20  used  vibrational  energy  of  deuterated  H20  reagents  to  overcome  the  large 
endothermicity  of  Reaction  (1), 

0(3P)  +  HOD(X,1Ai,4vOH)  — »  OH(X,2Ti)  +  OD(X,2Ti).  AH  =  +16.9  kcal  mol'1  (3) 

In  these  experiments,  product  internal  state  distributions  were  measured  directly  through  laser- 
induced  fluorescence,  and  therefore  the  usual  ambiguities  in  interpreting  the  overlapping  spectra 
of  Reactions  (1)  and  (2)  were  removed.  These  measurements  yield  rich  dynamical  information 
on  the  surfaces  involved  in  Reaction  (1). 

For  completeness,  we  note  that  collisions  of  0(  P)  +  H20(X,  Ai)  at  very  large  relative 
collision  velocities  can  produce  electronically  excited  OH(A,2E+), 

0(3P)  +  H20(X,1Ai)  ->  OH(A,2Z+)  +  OH(X,2n)  ,  AH  =  +1 10.5  kcal  mol'1  (4) 

which  has  been  observed  in  space-based11  and  laboratory  molecular  beam  measurements.10  The 
electronically  excited  OH  is  apparently  formed  through  high-lying  conical  intersections  in  the 
potential  energy  surfaces  of  the  four-atom  complex.  The  surfaces  involved  in  Reactions  (1-2) 
are  part  of  a  large  manifold  of  electronic  states  which  give  rise  to  the  products  of  Reaction  (4). 
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In  this  paper,  we  present  global  potential  energy  surfaces  for  the  three  lowest  triplet  adiabatic 

surfaces  involved  in  0(  P)  +  H20(X,  Ai)  collisions  and  present  results  for  classical  dynamical 

calculations  on  these  surfaces.  The  main  purpose  of  the  paper  is  to  investigate  the  role  of 

Reaction  (1)  in  producing  spectral  radiation  in  the  1-5  p  region.  By  comparing  to  available 

measurements,  we  will  assess  the  validity  of  these  surfaces  for  further  dynamical  studies 

including  Reaction  (2)  in  order  to  form  a  complete  picture  of  0(  P)  +  HoO(X,  Ai)  collisions  over 

a  broad  energy  range.  The  present  work  is  also  of  interest  as  a  guide  to  the  design  and  analysis  of 

possible  future  hyperthermal  O-atom  +  H2O  beam  measurements.  In  these  experiments,  it  should 

be  possible  to  measure  the  OH  products  and  their  internal  state  distributions  directly,  rather  than 

through  spectral  emission,  over  a  range  of  energies  near  64  kcal  mol'1  (~8  km  s'1  relative 
2 1 

collision  velocity).- 
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The  surfaces  we  report  for  0(  P)  +  H20(X,  Ai)  are  spline-based  fits  based  on  ~20,000  fixed 
geometry  ab  initio  calculations  at  the  CASSCF+MP2  level  with  a  0(4s3p2dlf)/H(3s2p)  one 
electron  basis  set.  The  surfaces  are  used  in  quasi-classical  trajectory  (QCT)  calculations  to 
compute  rate  constants  of  reaction  (1)  as  a  function  of  temperature.  The  computed  rate  constants 
compare  very  well  to  available  measurements.  We  also  compute  the  total,  ro-vibrationally 
resolved,  and  differential  angular  cross  sections  for  reaction  (1)  with  QCT  on  the  fitted  surfaces 
at  fixed  velocities  from  threshold  at  ~4  km  s'1  to  1 1  km  s'1  relative  collision  velocity  (122.5  kcal 
mol'1).  We  compare  the  computed  cross  sections  to  available  space-based  and  laboratory  data. 
We  also  compute  cross  sections  for  reaction  (3)  and  compare  these  results  to  results  of  Ref.  20 
using  the  fitted  surfaces. 

A  major  finding  of  the  present  work  is  that  the  ro-vibrationally  excited  OH(X  II)  product 

3  1 

from  0(  P)  +  H20(X,  Ai)  collisions  are  a  significant  and  perhaps  dominant  contributor  to  1-5  p 
spectral  emission  above  ~40  kcal  mol'1  collision  energy.  This  finding  is  in  contrast  to  most 
interpretations  of  the  experimental  data  which  attribute  this  radiation  to  collisional  excitation  of 
H20(X,  Ai)  exclusively.  Another  important  result  is  that  OH(X“IT)  products  of  reaction  (1)  are 
formed  in  two  distinct  ro-vibrational  distributions.  The  ‘active’  OH  products  are  formed  with  the 
reagent  O-atom,  and  their  ro-vibrational  distributions  are  extremely  hot.  The  remaining 
‘spectator’  OH  is  relatively  ro-vibrationally  cold.  Summed  over  both  OH  products,  below  ~50 
kcal  mol'1  collision  energy,  vibration  dominates  the  OH  internal  energy.  Above  ~50  kcal  mol"1, 
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rotation  is  greater  than  vibrational  energy,  but  the  total  product  energy  fraction  is  mostly  in 
translation.  We  note  that  the  present  surfaces  can  also  be  used  to  investigate  collisional  excitation 
ofH20(X,1Ai),  reaction  (2),  or  OH(X,2II)  +  OH(X,2n)  collisions,  the  reverse  of  reaction  (1),  but 
we  leave  such  work  for  a  separate  study. 

The  methods  used  for  the  potential  fit,  the  analytic  representation  of  the  fit,  the  fitted  surface 
properties  for  all  three  triplet  surfaces,  and  the  QCT  dynamics  methods  used  are  discussed  in 
Section  2.  Sections  3  presents  the  results  of  the  QCT  calculations  and  compares  results  to 
available  experimental  data.  Section  4  discusses  possible  future  work. 

II.  METHODS 

A.  Potential  Surfaces 

T.  I 

Collisions  of  0(  P)  with  H20(X  Ai)  give  rise  to  three  electronic  states  which  adiabatically 

2  2 

lead  to  OH(zn)  +  OH(  n).  These  three  lowest  triplet  adiabatic  states  are  the  focus  of  the  present 
study,  and  we  label  them  ‘state  1’,  ‘state  2’,  and  ‘state  3’,  where  state  1  is  the  lowest  energy  state 

and  state  3  is  the  highest.  We  fit  these  adiabatic  states  independently.  The  AH  of  the  0(  P)  + 

12  2  1 
H20(‘Ai)  ->  OHQIT)  +  OH(  IT)  reaction  is  +16.9  kcal  mol  .  From  rate  constant  measurements, 

the  activation  energy  is  — 17.0  kcal  mol"1  (Refs.  22,23)  or  a  relative  collision  velocity  of  4.1  km  s" 

,  which  implies  a  very  small  energy  barrier  of  ~0.1  kcal  mol"1.  We  note  that  crossings  of  the 

three  lowest  triplet  states  with  higher  triplet  states  can  lead  to  electronically  excited  OH(A“Z  ) 

with  an  endothermic ity  of  +  1 10.5  kcal  mol"1  (Ref.  12),  but  we  do  not  consider  such  interactions 

here.  Nor  do  we  consider  crossings  with  excited  singlet  oxygen,  O('D),  +  H20(X  1 A  i  )  states 

2  2  24 

which  can  also  lead  to  OH(“IT)  +  OH(  II).  We  further  note  that  by  symmetry  there  is  a  fourth 
triplet  leading  adiabatically  to  OH(  II)  +  OH(  II),  but  this  state  must  originate  from 
electronically  excited  triplets  of  O  +  H20.  Furthennore,  0(  P)  +  H20(X  AQ  reaction  channels 
leading  to  H2('Zg+)  +  02(3Zg"),  AH=-1.6  kcal  mol"1,  and  H(2S)  +  OOH(2A’),  AH  =+51.0  kcal  mol" 
',  are  not  addressed  here.  These  channels  go  through  high  lying  transition  states  (60-70  kcal  mol" 
')  7  and  should  be  accessible  at  hyperthermal  energies.  Fully  global  surfaces  for  triplet  O  +  H20 
should  take  these  channels  into  account.  A  summary  of  the  O  +  H20  reagent,  transition  states 
based  on  the  work  of  Ref.  17,  and  product  states  is  shown  in  Figure  1.  In  the  present  study  we 
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follow  the  three  lowest  triplet  states  of  0(3P)  +  H20(’Ai)  adiabatically  to  OH(2n)  +  OH(  fl) 
products. 

The  approach  taken  to  fit  the  three  lowest  triplet  adiabatic  surfaces  follows  closely  the  work 
of  Wu  et.  al.“  for  the  reaction  H  +  H2O  — >  OH  +  H2,  but  with  some  modifications  which  we 
discuss  below.  For  each  triplet  surface,  we  treat  the  full  (3N-6)  =  6  dimensions  by  defining  three 
‘active’  internal  coordinates  which  most  affect  motion  along  the  reaction  path  (bond  breaking 
and  forming)  and  three  remaining  ‘spectator’  internal  coordinates.  The  active  coordinate 
potential  is  a  3-dimensional  spline  detennined  from  a  large  set  of  ab-initio  calculations  spanning 
the  range  of  active  coordinate  values.  The  spectator  potential  is  determined  for  a  fixed  set  of 
active  coordinates  at  a  geometry  near  the  transition  state.  An  approximation  is  made  that  the 
spectator  coordinate  potential  does  not  vary  much  along  the  reaction  pathway  (changes  in  the 
active  coordinates).  A  global  analytic  representation  of  a  particular  surface  is  defined  as  a  sum  of 
active  and  spectator  parts.  The  global  potential  surface  connects  reagents  to  products  smoothly 
and  includes  the  proper  symmetrization  with  respect  to  O-atom  attack  on  either  of  the  H2O  H- 
atoms.  The  active  potential  spline  coefficients,  spectator  potential  parameters  and  other 
parameters  which  determine  the  surface  analytic  representation  are  based  on  new  ab-initio 
calculations  of  the  0(3P)  +  H20(!Ai)  potential  at  the  CASSCF+MP2  level.  The  internal 
coordinates,  analytic  potential  representation,  ab-initio  calculations,  and  properties  of  the  fitted 
surfaces  are  described  below. 

A.  1 .  Internal  Coordinates 

Figure  2  shows  the  internal  coordinates  used.  The  active  coordinates  are  the  set  {x,y,z} 
which  most  affect  reaction  going  to  OH  products.  The  coordinate  x  is  the  OH  bond  of  H2O  being 
broken,  y  is  the  OH  bond  being  formed,  and  z  is  the  angle  formed  from  between  the  breaking  and 
forming  bond,  atoms  Oa-Hb-Ob.  The  spectator  coordinates  are  the  set  {u,9,t},  where  u  is  the  OH 
bond  length  of  H20  which  is  not  attacked  by  the  incoming  O,  9  is  the  Ha-Oa-Hb  bond  angle,  and 
t  is  the  ‘out-of-plane’  torsion  angle  defined  as  the  angle  between  Ha  and  Ob  when  viewed  along 
the  Oa-Hb  bond.  The  active  coordinates  have  the  range  {x:  9— »oo,y:  9— »oo,z:  9— >189°},  while 
the  spectator  coordinates  have  the  range  {u:  9— >■  00 ,9:  9— >369°,  t:  9— >189°}. 
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A.  2.  Analytic  Representation  of  the  Surfaces 


We  first  define  the  “active+spectator”  potential,  Va+S  as  a  sum  of  active,  Va,  and  spectator,  Vs, 
parts, 


K+s  (*>  y, z > o, r)  =  ^  (*»  ^  (M>  0), 


Wo’Wo  ' 


(1) 


The  parameters,  (xo,yo,zo,Uo, 0o,To),  are  particular  values  of  the  internal  coordinates  defined 
further  below.  The  active  potential  is  further  broken  down  into  a  weighted  sum  of  two  parts  that 
depend  on  x, 

K  =  ^)-^^(x,y,z)uA.r^+^-W(T))-V^(x,y,z)llAir=m  .  (2) 

The  potential  Vspnne  is  defined  as  a  3D  spline  interpolation26  on  a  grid  of  calculated  ab-initio 
points  in  (x,y,z)  described  further  below.  Notice  that  the  spline  potential  only  depends  on  planar 
geometries  of  the  O  +  H2O  system  and  for  purposes  of  the  spline  interpolation,  z  is  considered  to 
be  continuous  from  0  to  360°.  W(z)  is  a  weighting  function  that  weights  the  x=0°  and  x  =  180° 
V ^1^  values  and  provides  an  interpolation  of  the  active  potential  for  non-planar  geometries, 
z  A  0°,180°  based  on  the  planar  Vspn„e  values.  W(r)  is  defined 

W(z)='£jamcosmz ,  (3) 

m— 0 


where  the  parameters  {am\  are  parameters  determined  from  a  fit  of  the  ab-initio  data  described 
below.  The  spectator  potential,  Vs,  is  a  sum  of  two  terms, 


v.=vm 


+VM*)\ 


(4) 


Here  Vu  is  a  Morse  potential  given  by, 
Vu(u)  =  D,[ 


(5) 


with  fit  parameters  {De,uo,b}.  The  other  part  of  the  spectator  potential,  Ve,  is  given  by, 
Ve  {0(x))  =  ( A  ■  (cos  6  -  cos  Gx  )4  +  B  ■  (cos  9  -  cos  9X  )3 )  ■  e ~a(x~Xi  1 
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(6) 


where 


0,  =  C  +  D  ■  e“a'(x“X|)  (7) 

with  fit  parameters  {A,B,C,D,a,a’,Xi}.  The  functional  form  of  Eq.  (6)  is  slightly  different  from 
the  study  of  Wu  et  al.  to  better  match  the  ab-initio  O  +  H20  surfaces. 

27  28 

Since  functional  forms  for  the  H20  reagents  and  OH  products"  have  already  been 
developed  which  match  experimental  data  very  well,  we  incorporate  these  known  potentials  in 
the  present  global  fit.  We  do  this  in  a  similar  fashion  to  Wu  et  al.  by  employing  switching 
functions  that  turn  on  and  off  the  reagent,  product,  and  Va+S  potentials  in  appropriate  regions.  We 


choose  the  following  forms  for  these  switching  functions: 

Sr  =  ~[1  +  tanh(r0  •  (y  -  q))]  (8a) 

Sq  =  ~  [1  +  tanh((y  -  x))]  (8b) 

SP  =|[l  +  tanh(/v(*-/h))]  (8c) 

with  fit  parameters  {r0,ri,p0,pi}.  Defining  Vr  to  be  the  isolated  H20  potential,  and  Vp  as  the 
potential  of  two  isolated  OH  products,  the  global  potential  can  then  be  expressed  as, 

=ST-S,-rT+^- Sr  •  S,  - Sp  •  (l - st ))•  F„,  +  Sp  ■  (l - S, )•  Vp  (9) 


An  undesired  side-effect  of  separating  active  and  internal  coordinate  potentials  is  that  the 
potential,  Vunsvmt  treats  the  hydrogen  atoms  on  the  water  reagent  in  an  unsymmetrical  fashion. 
This  means  that  the  potential  VunSym  is  different  depending  on  whether  the  incoming  oxygen 
attacks  Hb  rather  than  Ha,  which  is  unphysical.  To  compensate,  we  follow  Wu  et  al.  and 
introduce  the  switching  function  Ss,  which  makes  attack  on  either  hydrogen  equivalent. 

S,  =  \  I1  -  lanh(-S0  •  Ku„  -  RObHb  E  (  1  °) 

where  {so}  is  a  fit  parameter.  The  full  global  symmetrized  potential  is  then  the  weighted  sum  of 
two  unsymmetrized  parts, 
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Kym  =  Ss  •  vunsym (H.  -  Oa  -  Hb  +  Ob )  +  (1  -  S, ) ■  Vunsym (Hb  -  Oa  -  Ha  +  Ob ) 


(11) 


where  Vunsym(Ua  -Oa  -Hb  +  Ob)  is  the  potential  of  the  atoms  in  their  original  configuration,  and 
+  Ob)  is  the  potential  with  the  H  atoms  switched.  This  is  essentially  an 
averaging  of  two  geometric  configurations. 

Table  I  gives  the  values  of  all  the  parameters  for  the  potential  for  all  three  electronic  states 
and  how  they  were  derived.  Many  of  the  values  were  obtained  from  ab-initio  calculations 
described  in  the  next  section. 


TABLE  I.  Surface  fit  parameters.  All  values  are  a.u. 


Fit  Parameters 

Values  (a.u.) 

How  Determined 

“Transition  State” 
{xo,yo,zo,uo,0oAo} 

(2.3501,2.0558,  141.06°, 1.8287,  104.55°, 0.0) 

Ref.  17 

W(t),  {a0,...,a6} 
state  1 
state2 
state3 

(  0.672,0.649,-0.311,-0.347,0.373,0.198,  -0.234) 
(1.47,-0.296,-0.387,0.511,-0.0724,0.263,-0.496) 
(3.63,  0.480,  -5.74,  -0.614,  5.15,  0.627,  -2.55) 

Fit  of  ab-initio  data 

Vu  {De,b,u0} 

state  1 

state2 

state3 

(0.1628,  1.2200,  1.8500) 

(0.1757,  1.2644,  1.8287) 

(0.1574,  1.3000,  1.8286) 

Fit  of  ab-initio  data 

Ve  {A,B,C,D,a,a’,xi} 

state  1 

state2 

state3 

(0.7400,0.065,99.0,  110.,  1.9,  1.6,2.11) 

(0.221,  0.0164,  98.7,  3.3,  1.0525,  0.8,  2.232) 

(0.137,  -0.0416,  103.5,  16.,  2.307,  0.1,  2.35) 

Fit  of  ab-initio  data 

switching  functions 
{r0,ri,po,piSo} 

(4.0,  6.0,  3.0, 4.0, 3.0) 

Ref.  25 

OH  product 

{De,b,u0} 

(0.1698,  1.21,  1.85) 

Ref.  28 

H20  reagent 

many 

Ref.  27 

V spline  coefficients 

many 

Fit  of  ab-initio  data 
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A.  3.  Electronic  Structure  Calculations 

Most  of  our  ab-initio  calculations  were  based  on  a  10  electron  8  orbital  (10e8o)  CASSCF 
+  MP2  (Complete-Active-Space  Self-Consistent-Field  +  Second  Order  Perturbation  Theory) 
treatment  with  a  (4s3p2dlf)  O-atom  basis  and  3s2p  H-atom  basis  for  a  total  of  90  basis 
functions.  The  (10e8o)  CASSCF  includes  all  excitations  of  eight  electrons  within  the  space  of  an 
initial  configuration  of  two  bonding,  two  lone-pair,  two  singly  occupied,  and  two  un-occupied 
anti-bonding  orbitals.  (The  Is  and  2s  oxygen  orbitals  are  not  included  in  the  CASSCF  space.) 
This  totals  1512  configuration  state  functions.  We  state-average  the  three  lowest  electronic  states 
and  perform  all  these  calculations  with  no  symmetry.  This  choice  of  CASSCF  space  ensures  that 
the  three  electronic  states  treated  will  be  degenerate  on  the  reagent  and  product  sides  of  the 
reaction  as  they  should  be,  and  that  the  CASSCF  space  includes  bonding  and  anti-bonding 
orbitals  so  that  OFI  bonds  can  break  and  form  in  an  even  handed  way  along  the  reaction  path.  A 
summary  of  the  electronic  energies  and  energy  differences  of  the  reagents,  products,  and  at  the 
lowest  transition  state  geometry  (as  determined  from  the  study  of  Ref.  17)  are  given  in  Table  II. 

9Q 

All  calculations  were  done  with  the  GAMESS  electronic  structure  package. 


TABLE  II.  Computed  CASSCF-MP2  Energies.  Te  is  defined  to  be  the  calculated  energy 
difference  with  respect  to  O  +  FEO(Re),  Energy-Energy(0+H20(Re)). 


Geometry 

Energy 

(a.u.) 

T 

1  e 

(kcal  mol'1) 

0  +  H20(Re) 

-151.289431 

0. 

OH(Re)  +  OH(Re) 

-151.264587 

15.59 

Transition  State 

(Following  Ref.  17,  Table  I) 

state  1 

-151.266063 

14.66 

state  2 

-151.250827 

24.22 

state  3 

-151.227117 

39.10 
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The  experimental  Te  leading  to  OH  +  OH  is  14.7  kcal  mol'1  which  is  0.9  kcal  mol'1  below  the 
present  calculated  ab  initio  results.  Because  we  use  switching  functions  to  smoothly  connect  to 
the  reagent  and  product  parts  of  the  potential,  we  have  defined  the  Te  for  the  present  fitted 
potential  to  be  the  experimentally  measured  Te  of  14.7  kcal  mol1.  Computing  transition  state 
geometries  at  the  CASSCF-MP2  level  is  problematic,  especially  with  state  averaging.  So  we 
have  adopted  the  transition  state  geometry  reported  in  Ref.  17  which  uses  a  basis  set 
extrapolation  technique  within  a  QCISD(T)  treatment.  Ref.  17  obtained  a  Te  for  the  first 
transition  state  of  state  1  of  +  1.61  kcal  mol'1  above  the  OH  +  OH  limit,  while  the  present  surface 
would  be  nearly  equal  in  energy  to  the  OH  +  OH  limit  for  the  lowest  triplet.  Deyerl  et  al.19  obtain 
a  lowest  triplet  transition  state  very  similar  in  geometry  to  Karkach  et  al.  within  a  QCISD(T)/6- 
31 1++G(3df,2p)  treatment  which  is  +0.10  kcal  mol’1  above  the  OH  +  OH  limit  (electronic 
energy  -151.315271).  As  discussed  further  below,  the  second  and  third  lowest  triplets  have 
potential  surfaces  quite  similar  to  the  lowest  triplet,  except  that  the  barrier  heights  are  larger  by 
-10  and  ~15  kcal  mol'1  respectively.  The  second  and  third  states  should  therefore  play  a  role  at 
hyperthermal  energies. 

To  obtain  the  spline  fit  parameters,  GAMESS  CASSCF+MP2  calculations  were  done  on  a 
34x34x14  grid  of  x,y,  and  z  points  for  a  total  of  16,184  points.  For  these  calculations  (0=104.5°, 
u=  1.8085  a.u.,  t=0°).  The  values  of  u  and  0  are  slightly  different  from  the  transition  state 
geometry  adopted  here,  but  these  differences  will  have  minimal  impact  on  the  fit.  The  x  and  y 
coordinate  grid  spanned  the  range  1.3228  a.u.  to  7.559  a.u.  with  spacings  of  0.1890  a.u.  in  each 
direction.  The  z-coordinate  values  extend  from  50°  to  310°  with  a  spacing  of  20°.  We  fit  the 
computed  energies  to  a  3D  cubic  spline  (with  tension)  after  Ref.  26.  For  geometries  lying 
outside  the  x-  and  y-value  grid  points,  we  map  the  results  to  the  nearest  extreme  grid  value.  For 
z-values  outside  the  calculated  range,  we  perform  a  linear  interpolation  between  the  nearest 
extreme  grid  values.  Also,  for  z  e  [50°,70°]  we  found  it  necessary  to  use  a  multi-dimensional 
linear  (rather  spline  fit)  to  avoid  un-physical  oscillations  in  the  interpolated  potential.  We  note 
the  cubic  spline  interpolation  is  simultaneously  performed  for  all  three  active  coordinates,  so  the 
z  coordinate  is  treated  on  equal  footing  with  the  x  and  y  coordinates.  We  also  note  that  for 
subsequent  dynamics  calculations,  splines  are  especially  efficient  for  evaluations  of  the  potential 
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as  all  interpolations  are  perfonned  locally  while  the  value  of  the  function  and  its  first  derivative 
remain  continuous  from  one  grid  point  to  the  next. 

To  obtain  the  W(r)  fit  parameters  for  each  state,  we  performed  CASSCF+MP2  calculations 
as  described  above  fixing  all  internal  coordinates  except  x  to  the  transition  state  geometry.  We 
then  varied  x  in  steps  varying  from  5  to  20  degrees  from  0  to  180.  Similarly,  to  obtain  the  Vu  fit 
we  performed  CASSCF+MP2  calculations  fixing  all  internal  coordinates  except  u  to  the 
transition  state  geometry  and  varied  u  over  61  points  in  the  range  1.22  to  9.44  a.u.  For  the  Ve  fit 
we  performed  a  series  of  CASSCF  +  MP2  calculations  fixing  all  internal  coordinates  except  9 
and  x  at  their  transition  state  geometries.  We  varied  9  and  x  on  a  two-dimensional  grid,  varying  9 
every  5  degrees  from  0  to  360  degrees  and  varying  x  from  1.16  to  8.019  a.u.  over  24  unevenly 
spaced  steps.  Since  Vspune  already  contains  a  dependence  on  x,  we  fit  the  residual  dependence, 
V ab-initio  ~  V spune  in  order  to  isolate  the  contribution  of  varying  the  spectator  coordinate  9. 

A.  4.  Properties  of  Potential  Surfaces 

The  potential  surface  fits  for  the  three  lowest  triplet  states  of  the  O  +  H2O  system  are  shown 
in  Figure  3  for  planar  geometries.  In  figure  3,  the  spectator  coordinates,  {u,9,x},  are  fixed  at  their 
transition  state  values  as  given  in  Table  I  to  isolate  the  Va  portion  of  the  potential,  and  the  z- 
coordinate  is  fixed  at  89°,  149°  and  249°.  The  active  coordinate  {x}  varies  along  the  abscissa  the 
active  coordinate  {y}  varies  along  the  ordinate.  Large  values  of  y  correspond  to  the  separated  O 
+  H2O  reagents  and  large  values  of  x  correspond  to  separated  OH  +  OH  products.  The  transition 
state  barrier  for  this  fitted  potential  occurs  near  z=149°  for  all  three  surfaces  and  its  energies  are 
— 15,  ~24,  and  ~39  kcal  mol,  respectively  with  respect  to  O  +  H20(Re).  The  transition  state 
saddle  point  for  the  lowest  state  occurs  near  x  «  2.4  and  y  ~  2. 1  a.u.,  so  that  the  OH  bond  that  is 
breaking  along  x  is  larger  than  the  OH  bond  that  is  forming  along  y. 

Figure  4  shows  the  fitted  potential  Va  for  non-planar  geometries  for  all  three  fitted  surfaces. 
In  Figure  4,  all  coordinates  are  fixed  at  their  transition  state  values  as  shown  in  Table  I,  and  the 
torsion  angle  x  is  varied  from  9  to  189°  to  isolate  the  effect  of  out-of-plane  motion.  For  the 
lowest  electronic  state  there  is  a  preference  for  9°  (cis)  geometries  over  189°  (gauche)  by  ~7  kcal 
mol"1.  There  appears  to  be  an  avoided  crossing  near  x=139°  between  the  two  lowest  electronic 


12 

Approved  for  public  release;  distribution  unlimited. 


states.  For  the  third  electronic  state  there  is  an  energy  maximum  near  t=90°  and  the  0  and  180° 
geometries  have  about  the  same  energy.  The  barriers  for  out-of-plane  motion  are  generally  small. 

Figure  5  shows  the  Ve  component  of  the  spectator  potential,  which  is  sensitive  to  the  reagent 
FFO  internal  bond  angle,  for  all  three  fitted  surfaces.  For  these  plots,  all  internal  coordinates 
except  9  are  fixed  at  their  transition  state  values  as  given  in  Table  I.  For  all  three  surfaces  there  is 
a  barrier  at  linear  H2O  geometries,  which  is  -30  kcal  mol'1  for  the  lowest  electronic  state.  There 
are  local  minima  at  9  «  194.5°  and  9  «  255.5°  which  correspond  to  cis  and  trans  geometries  of 
the  combined  O  +  H2O  system  where  the  reagent  is  close  to  its  isolated  equilibrium  bond  angle. 
Ve  is  constructed  to  be  zero  here  as  the  Va  part  of  the  potential  already  treats  these  geometries.  At 
small  and  large  values  of  9  correspond  to  ‘small  angle’  geometries  of  the  reagent  FFO,  and  so  are 
dominated  by  the  nuclear-nuclear  repulsions  of  the  FFO  hydrogens.  Figure  5  is  for  fixed  a  value 
of  x.  As  x  increases  and  products  are  formed,  Ve  retains  the  same  shape  but  decreases  in 
magnitude  for  all  three  states  as  expected. 

Figure  6  shows  the  full  potential  for  all  three  fitted  states,  Vsym  for  planar  geometries.  In 
Figure  6,  the  coordinates  {x,u,9,x}  have  been  fixed  at  their  transition  state  values  as  shown  in 
Table  I,  and  {y,z}  have  been  allowed  to  vary.  This  view  shows  the  potential  for  a  fixed  H20 
reagent  at  varying  angles  of  approach  of  the  O-atom.  The  chief  feature  of  these  surfaces  is  a  well 
near  angles  of  approach  corresponding  to  cis  geometries  of  the  O  +  H2O  system,  in  regions 
centered  near  z  «149°  for  the  lowest  electronic  state.  For  the  higher  electronic  states  the  well  is 
not  so  deep  and  z  «189°.  These  z  angles  are  the  preferred  angles  of  approach  of  the  O-atoms  for 
chemical  reaction.  As  the  reagent  OH  bond  lengthens  (x  increases)  and  the  product  OH  bond 
forms  (y-decreases)  the  shallow  wells  in  Figure  6  deepen.  Generally,  other  angles  of  approach 
are  repulsive,  except  at  large  distances  where  there  are  shallow  long-range  type  interactions  of 
the  approaching  Oxygen  and  H2O  hydrogens.  Some  of  this  long-range  attraction  corresponds  to 
geometries  affected  by  our  separation  of  active  and  spectator  coordinates  and  the  way  the  present 
potential  is  symmetrized,  and  so  we  do  not  expect  these  regions  to  be  as  accurate  as  the  transition 
state  geometry  regions  which  are  fit  by  splines  alone. 


II.  B.  Quasi-classical  trajectory  (QCT)  calculations 
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31* 

We  use  the  VENUS  classical  trajectory  code  with  an  added  interface  modification  in  order 
to  use  our  fitted  potential  surfaces  and  derivatives  to  generate  rate  constants  and  cross  sections. 
For  all  calculations  we  use  a  fixed  time  step  of  8t  =  l.e-16  seconds  and  a  maximum  impact 
parameter  bmax  ranging  from  5.7  a.u.  to  7.56  a.u.  To  generate  rate  constants  at  fixed  temperatures, 
we  perform  l.e5  to  5.e5  trajectories  per  electronic  state  per  temperature.  To  generate  cross 
sections  at  fixed  translational  energies,  we  fix  the  reagent  water  in  its  ground  ro-vibrational  state 
and  perform  l.e5  to  5.e5  trajectories  per  each  electronic  state  per  energy.  To  obtain  total  rate 
constants  or  cross  sections  we  add  the  trajectory  contributions  from  each  electronic  state  together 
and  multiply  by  1/3,  the  electronic  degeneracy  factor.  For  final  state-resolved  quantities  the 
OFI(v,j)  states  are  binned  in  the  usual  way  as  outlined  in  Ref.  3 1 .  Final  FEO  products  are  binned 
according  to  the  amount  of  internal  energy  in  rotation  and  vibration.  We  note  that  our  results  do 
not  adjust  for  zero-point  energy  of  the  reagents,  which  will  sometimes  tend  to  make  such 
classical  dynamical  results  more  ‘reactive’  than  the  potential  actually  allows. 

III.  RESULTS  AND  DISCUSSION 

A.  Rate  constants  at  fixed  temperatures  for  0(3P)  +  H20(X,1Ai)  — >  OH(X,2IT)  +  OH(X,2II) 

3  1 

Figure  7  shows  the  computed  QCT  rate  constant,  krxn,  for  the  reaction,  0(  P)  +  H20(X,lAi) 

2  2 

->  OH(x/n)  +  OH(X,  FI),  as  a  function  of  temperature  together  with  data  from  two 

22  23 

measurements.  ’  Calculations  were  performed  from  1,000  K  to  3,500  K  at  intervals  of  500  K. 
There  is  generally  good  agreement  between  the  measurements.  (especially  Ref.  23), 
Generally,  the  present  results  are  within  experimental  uncertainties  (about  a  factor  2).  Usually, 
neglect  of  a  zero  point  energy  correction  as  in  the  present  QCT  calculations  would  tend  to  over¬ 
estimate  the  rate  constant,  but  there  may  be  a  compensating  error  in  the  neglect  of  a  tunneling 
correction.  We  do  not  report  results  below  1,000  K,  because  of  the  difficulty  of  obtaining 
acceptable  QCT  statistics.  The  separate  electronic  state  contributions  to  the  QCT  computed  rate 
constant  are  also  shown.  The  dominant  contribution  is  from  the  lowest  energy  electronic  state, 
state  1 ,  as  expected.  The  second  lowest  electronic  state  also  makes  an  appreciable  contribution  at 
all  temperatures  studied  and  at  the  highest  temperatures  has  a  contribution  within  a  factor  of  two 
of  state  1 . 
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B.  Cross  sections  at  fixed  collision  energies  for  0(3P)  +  H20(X,1Aj)  -»  OH(X,2II)  + 

oH(x,2n> 

1.  Total  and  vibrationally  resolved  cross  sections 

Figure  8  shows  the  total  calculated  QCT  ‘OH  cross  section’,  <JOH,  for  the  reaction,  0(  P)  + 

1  2  2 

H20(X,  Ai)  ->  OH(X,zn)  +  OH(X,  FI),  as  a  function  of  collision  energy  along  with  component 
contributions  from  the  three  lowest  triplet  electronic  states.  Calculations  were  perfonned  at  the 
following  collisions  energies:  16.2,  18.7,  20.5,  25.3,  30.6,  36.4,  49.6,  64.8,  82.0,  101.2,  and 
122.5  kcal  mol'1  which  correspond  to  relative  collision  velocities  of  4.0,  4.3,  4.5,  5.0,  5.5,  6.0, 
7.0,  8.0,  9.0,  10.0,  and  1 1  km  s'1.  We  define  the  OH  cross  section,  aOH ,  in  the  following  way, 

&OH  =  2(7rXn  =  ^max  ^  >  (  1 2) 

where  urxn  is  the  reaction  cross  section,  Nr  is  the  number  of  reactive  QCT  trajectories  and  Ntot  is 
the  total  number  of  trajectories.  The  cross  section,  aOH ,  is  more  convenient  for  comparisons  to 

most  of  the  available  measurements  which  often  report  absolute  cross  sections  by  detecting  and 
adding  together  signals  from  both  OH  products  formed  in  the  reaction.  The  OH  cross  section  is 
relatively  small  near  threshold  and  increases  with  energy.  By  30  kcal  mol'1  collision  energy,  the 
cross  section  is  relatively  constant,  although  it  continues  to  increase  with  collision  energy.  This 
implies  that  future  measurements  need  to  start  below  5  km  s'1  to  capture  the  rapid  increase  of  the 
excitation  function  near  threshold.  As  the  energy  increases  past  the  energy  barriers  of  state  2  and 
state  3,  they  begin  to  contribute  significantly  to  the  cross  section.  At  the  highest  collision 
energies  studied,  the  contribution  from  state  3  becomes  comparable  to  that  of  state  2,  where  they 
are  about  2.5  times  smaller  than  state  1 . 

Figure  9a  shows  the  QCT  vibrationally  resolved  OH  cross  sections,  oOH{y),  for  the  reaction, 

0(3P)  +  H20(X,1A1)  ->  OH(X,2n)  +  OH(X,2n),  at  25.3,  49.6,  64.8,  82.0,  and  122.5  kcal  mol'1 
collision  energy  or  5.0,  7.0,  8.0,  9.0,  and  11.0  km  s'1  collision  velocity,  respectively,  where  we 
have  summed  the  contributions  from  both  OH  products  in  a  particular  vibrational  state  (Eq.  12). 
At  5,7,  and  8  km  s'1  collision  velocity,  the  distributions  are  well  fit  by  thennal  populations  with 
temperatures  of  1,900  K,  2,700  K,  and  3,100  K  respectively.  At  9  and  11  km  s'1  collision 
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velocity  we  can  fit  the  lower  OH(v)  states  both  with  thermal  temperatures  of  -3,500  K.  The 
highest  OH(v)  populations  for  these  higher  velocities  begin  to  show  a  different  type  of 
distribution.  However  these  cross  sections  are  very  low,  and  proper  fitting  of  these  population 
requires  better  statistics.  Figure  9b  and  9c  show  the  component  ‘active’  and  ‘spectator’ 
vibrationally  resolved  OH  cross  sections.  The  active  OH  products  are  from  the  incoming  reagent 
atomic  oxygen,  and  the  spectator  OH  products  are  made  from  the  oxygen  in  the  reagent  H2O,  so 
that  the  total  OH  cross  section  (Fig.  9a)  for  a  given  vibrational  level  is  the  sum  of  active  and 
spectator  OH  cross  sections  for  that  level.  Although  measurements  of  the  active  and  spectator 
OH  cross  sections  would  probably  only  be  available  through  detailed  molecular  beam  isotopic 
studies,  the  computed  active  and  spectator  OH  cross  sections  offer  insight  into  the  dynamics  of 
the  collision  process.  For  example,  it  is  clear  that  the  active  and  spectator  vibrationally  resolved 
OH  cross  sections  are  quite  distinct.  The  spectator  cross  sections  are  significantly  colder  than  the 
active  OH  populations.  The  spectator  OH’s  mainly  populate  OH(v’=0,l)  levels,  with  the  v’=l 
relative  populations  significantly  less  than  that  of  the  active  OH(v’=l)  populations.  The  spectator 
OH(v’>l)  populations  are  almost  negligible.  This  means  that  the  total  OH(v’>0)  populations  are 
due  almost  exclusively  to  the  active  OH,  and  the  total  OH(v=0)  populations  are  mainly  due  to  the 
spectator  OH.  The  OH  bond  being  fonned  in  the  reaction  takes  up  most  of  the  vibrational  energy 
in  the  products.  Joint  active/spectator  distributions  would  reveal  more  detailed  information,  but 
would  require  substantially  better  statistics. 

The  vibrationally  resolved  results  of  Figure  9a  can  be  compared  with  data.  Fixed  energy 
spaced-based  and  laboratory  measurements  for  0(3P)  +  H20(X,1Ai)  — »  OH(X,2n)  +  OH(X,2n) 
focus  on  the  detection  of  infrared  radiation  from  vibrationally  excited  total  OH  (active  + 
spectator,  Fig.  9a)  products  summed  over  the  populated  vibrational  states.  Since  the  available 
spaced-based  and  laboratory  measurements  do  not  measure  vibrationally  resolved  cross  sections 
directly,  it  is  necessary  to  briefly  review  the  nature  of  these  experiments  to  make  appropriate 
comparisons  with  the  present  results. 

As  outlined  in  Ref.  7,  space-based  measurements  (such  as  those  of  Refs.  5,  6  and  9)  are 
characterized  by  a  steady-state  outgassing  or  exhausting  source  of  H20  with  a  known  or  assumed 
rate,  NHO  (number  of  H20  per  second),  interacting  with  atomic  oxygen  in  the  atmosphere.  The 

signal  is  usually  observed  through  a  somewhat  opaque  atmosphere  integrated  over  a  wavelength 
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band  and  characterized  by  a  transmittance,  Tax-  The  observed  steady  state  intensity,  Is^ce 
(photons/second),  can  be  written, 


j  space 
^  AA 


G 


G1r 


N  T 

iV  H.O1  AA  9 


(13) 


where  om  is  the  total  (reactive  and  non-reactive)  collision  cross  section  for  O  +  H20  and  a*  is 
the  photon  excitation  cross  section  which  includes  all  contributions  due  to  photon  cascades  from 
multi-quantum  excited  states.  For  detection  in  the  2-5  p  spectral  region  where  most  of  the 
relevant  optical  measurements  take  place  for  this  system,  the  intensity  will  be  sensitive  to  single 
vibrational  quantum  changes,  so  we  write  a*  as, 


(i4) 

species  v=l 

where  v  is  the  vibrational  quantum  number,  (j'pecies  is  the  cross  section  to  a  particular  vibrational 

state  of  either  collisionally  excited  Ff20  or  reaction  product  OH  radiating  in  the  detection  band.  It 
is  important  emphasize  that  observations  usually  include  contributions  from  both  species,  and  so 
the  first  sum  extends  over  OH  and  H20.  The  second  sum  extends  over  all  possible  excited 
vibrational  levels  of  a  given  species.  Space-based  measurements  usually  assume  or  otherwise 
employ  models  to  estimate  values  for  N HO , 7)v  and  atot,  and  from  the  observed  intensity,  V^ce , 

report  a  photon  excitation  cross  section,  a*.  Therefore,  uncertainties  in  the  estimated  values  of 
N H  O ,  Ta.  and  atot  are  reflected  in  the  reported  a*. 

On  the  other  hand,  ground-based  laboratory  measurements  involving  optical  detection 
typically  use  pulsed  beam  sources.  Signals,  I1^  ,  from  such  measurements  are  related  to  the 
photon  excitation  cross  section  in  the  following  way, 

iTx «  Z  YaATcies <pecie\  (is) 

species  v=l 

where  A°pecws  are  Einstein  coefficients  (photons/second)  for  single  photon  radiative  decay. 

Because  contributions  from  separate  species  and  from  individual  vibrational  levels  usually 

cannot  be  distinguished,  lab  measurements  typically  report  a  single  value  for  the  right  hand  side 
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double  sum  of  Eq.  (15)  at  fixed  O  +  H2O  collision  velocities,  and  for  convenience  we  denote  this 
quantity  as  Act (lengtir/second).  We  note  that  these  laboratory  measurements  are  often  difficult 
to  carry  out  and  calibrate  absolutely  with  high  precision  due  to  uncertainties  in  reactant 
concentrations  and  the  inability  of  optical  systems  to  capture  all  collision  produced  photons 
within  a  limited  field  of  view. 

Table  III  gives  two  spaced-based  results  for  the  photon  excitation  cross  section  (Eq.  14),  two 
laboratory  results  for  the  product  cross  sections  times  Einstein  A  (Eq.  15),  and  the  corresponding 
quantities  derived  from  the  present  calculated  vibrationally  resolved  total  OH  cross  section  (Fig 
9a).  These  quantities  follow  the  analysis  of  Ref.  7  and  assume  A°H  =  v  •  17  sec  1 .  It  is  estimated 

that  the  experimental  results  are  uncertain  by  a  factor  of  ~2  to  ~  3.  At  4.3  km  s'1,  the  laboratory 
results  of  Ref.  1  are  about  12  times  larger  than  the  presently  computed  contribution  from  excited 
OH.  This  collision  energy  is  too  low  to  excite  OH  populations.  Therefore,  the  experimental 
signal  is  nearly  all  H2O,  collisionally  excited  in  the  Vi  and  V3  stretching  modes.  At  8  km  s'1,  well 
above  the  reaction  energy  barrier  to  make  OH,  the  measured  space-based6  results  are  about  a 
factor  of  2  smaller  and  laboratory  results4  about  equal  to  the  presently  computed  OH  contribution 
to  the  overall  signal.  At  1 1  km  s'1  there  is  good  agreement  between  the  space  based  combined 
signal  and  presently  computed  contribution  from  OH.  The  present  results  therefore  imply  that  at 
and  above  8  km  s'1  collision  velocity,  the  OH  contribution  is  at  least  comparable  and  perhaps 
larger  than  the  H20  contribution  to  the  observed  infrared  signals.  We  note  that  the  H2O 
collisional  excitation  cross  sections  have  been  computed  in  previous  studies.  Although  these 
studies  used  the  same  non-reactive  potential  surface,  the  quantum  scattering  results  of  Ref.  14 
differ  by  orders  of  magnitude  from  the  classical  scattering  results  of  Ref.  13  even  at  high 
velocities,  so  it  is  difficult  to  assess  the  possible  relative  contributions  from  collisionally  excited 
water.  We  note  however,  that  at  8  km  s'1  collision  velocity  there  is  fairly  good  agreement 
between  the  quantum  results  for  collisional  excitation  of  H2O  and  the  measurements  of  Ref.  4. 
With  the  presently  computed  contribution  from  OH,  the  theoretically  derived  combined  H2O/OH 
would  still  be  within  the  experimental  error  bounds. 
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TABLE  III.  Excitation  cross  sections  derived  from  spaced-based,  laboratory,  and  presently 

3  1 

calculated  results  for  0(  P)  +  H20(X  Ai)  collisions  to  produce  radiation  from  2-5  p. 


Relative  Collision 

velocity  (km  s'1) 

V=1 

X  Z^r" 

species  v=l 

ZLA” 

V=1 

species  v=l 

(Theory:  a.u.) 

(Exp:  a.u.) 

(Theory: 

a.u./second) 

(Exp.  a.u./second) 

4.3 

0.09  (Ref.  1) 

1.07 

~8. 

1.28 

0.54  (Ref.  6) 

21.75  (Ref.  4) 

18.20 

11 

4.41 

3.57  (Ref.  9) 

We  note  that  the  vibrational  temperature  of  3,100  K  obtained  from  the  present  calculations  at 
8.0  km  s"1  is  quite  different  from  what  has  been  conjectured  by  Bernstein  et  al.7  in  their  fit  of 
observed  spectra  of  Ref.  4  obtained  in  hyperthermal  O-atom  beam  experiments  at  similar 
collision  energies.  Bernstein  et  al.7  derive  a  vibrational  temperature  of  7,000  K,  which  was 
mostly  driven  by  the  observed  large  OH  overtone  signals.  The  present  results  cannot  reproduce 
the  apparently  large  overtone  populations  observed  in  Ref.  4. 


2.  Rotationally  resolved  cross  sections 

Figure  10  shows  the  rotationally  resolved  total  OH  cross  at  5.0,  7.0,  8.0,  9.0,  and  1 1.0  km  s'1 
relative  collision  velocity  for  0(3P)  +  H20(X,'Ai)  — »  OH(X,2n)  +  OH(X,2n),  where  the  OH(v’) 
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populations  have  been  summed  over.  These  cross  sections  show  a  bi-model  distribution,  a 
narrow  low  OH(j’<10)  peak  and  a  wide  almost  flat  OH(j’)  extending  to  very  high  OH(j’)  which 
eventually  falls  off  as  the  energetic  limit  is  reached.  As  shown  in  Fig.  11  for  8  km  s'1  collision 
velocity,  the  two  underlying  rotational  distributions  come  from  the  active  and  spectator  OH’s. 
The  active  OFI  is  much  more  rotationally  excited  than  the  spectator  OFI,  so  that  the  active  OFI 
receives  the  vast  majority  of  the  internal  energy  imparted  to  products  in  the  collision.  The 
superimposed  oscillatory  structure  for  low  active  OH(j)  populations  appears  to  be  real  and  not 
caused  by  statistical  uncertainty. 

We  have  fit  the  separate  active  and  spectator  OH(j  ’)  rotational  populations  to  a  distribution 

32 

derived  from  surprisal  theory,  ' 

P(vj)  oc  exp {6rEvj  )( 2j  +  1)2(2/  +  1  )(E  -  Evj  -  Eml  )‘/2 ,  (16) 

ml 

where  (y /)  are  the  vibrational  and  rotational  quantum  numbers,  respectively,  for  OH  product 
being  fit,  ( ml)  are  the  vibrational  and  rotational  quantum  numbers,  respectively,  for  the  other  OH 
product,  E  is  the  total  energy  available  to  be  distributed  to  products  in  the  reaction,  6,  is  a 
surprisal  parameter,  the  sum  of  {ml)  extends  over  all  the  available  energy  states,  and  we  have  set 
v=0.  The  results  are  shown  in  Fig.  1 1  for  separate  fits  to  the  active  and  spectator  OH’s,  where  6r 
=  0.0  (statistical)  for  the  active  OH  and  0r  =  -0.41  for  the  spectator  OH.  We  note  that  the  peak 
OH(j)  products  correspond  to  rotational  temperatures  of  1 1,400  K  and  1,100  K  for  the  active  and 
spectator  OH’s  respectively.  The  analysis  of  Ref.  7  of  the  data  in  Ref.  4  derives  a  rotational 
temperature  for  total  OH  of  1,500  K.  This  temperature  should  be  dominated  by  the  active 
vibrationally  excited  OH  for  which  we  calculate  a  much  higher  rotational  temperature. 


3.  Product  energy  disposal 

Figure  12a  shows  the  total  OH  and  Figure  12b  shows  the  active  and  spectator  OH  energy 
disposal  shown  as  energy  fractions  of  translation,  vibration,  and  rotation,  as  a  function  of  the 
collision  energy.  We  define  the  total  energy,  E  =  Econ  +  Ez.p.(H20)  -  14.7  kcal  mol'1,  and  the 
energy  fractoins  as  fv  =  <EV>/E,  fr  =  <Er>/E,  and  fj  =  l-<fv>-<fi>.  We  note  that  the  average 
energies,  <Ev>  and  <Er>,  were  obtained  before  histogram  binning  and  are  summed  over  all 
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trajectories  with  the  appropriate  weighting  of  electronic  state  contributions.  For  the  total  OH 
energy  disposal  at  energies  below  —35  kcal  mol'1  the  vibrational  and  translation  fractions  are 
about  the  same  and  are  appreciably  larger  than  the  rotational  fraction.  However,  above  50  kcal 
mol"1  the  rotational  fraction  begins  to  dominate  over  the  vibrational  fraction  as  the  vibrational 
energy  becomes  increasingly  diverted  to  translational  energy  with  increasing  collision  energy. 
Figure  12b  shows  the  breakdown  of  the  energy  fractions  into  active  and  spectator  components.  It 
is  quite  clear  that  the  active  OH  has  the  majority  of  product  internal  energy  and  this  tendency 
increases  dramatically  above  —40  kcal  mol1  collision  energy.  The  largest  changes  with  collision 
energy  occur  with  the  active  and  spectator  vibrational  energy  fractions  which  decrease 
dramatically  above  25  kcal  mol"1  collision  energy.  As  the  collision  energy  increases,  the 
vibrational  energy  fraction  becomes  diverted  into  rotation  and  translation,  especially  at  the 
highest  collision  velocities. 

4.  Differential  angular  cross  sections 

The  differential  angular  cross  section  summed  over  all  OH  states  is  shown  in  Fig.  13  at  5.0, 
7.0,  8.0,  9.0,  and  11.0  km  s"1  collision  velocity.  At  low  collision  velocities,  the  distribution  is 
backward  peaked,  but  as  the  collision  velocity  increases  the  angular  distribution  becomes  more 
and  more  forward  peaked,  until  at  1 1  km  s"1  products  are  mostly  forward  scattered.  There 
appears  to  be  a  change  in  mechanism  from  backward  to  forward  peaked  with  collision  velocity. 
This  mirrors  the  changes  in  energy  deposition  from  vibration  to  rotation  and  translation  with 
collision  energy,  and  indicates  different  reaction  mechanisms  at  low  and  high  energies.  The 
changeover  occurs  between  5  and  7  km  s"1. 


C.  Comparisons  to  measurements  of  0(3P)  +  HOD(X,4vOH)  — »  OH(X,2IT)  +  OD(X,2II) 

The  same  potential  surfaces  have  been  used  to  examine  the  measurements  of  Ref.  20  of  the 
reaction  0(3P)  +  HOD(X,4vOH)  — »  OH(X,2n)  +  OD(X,2n).  In  these  measurements,  atomic 
oxygen  formed  from  photolysis  of  NO2  collides  at  collision  energies  of  1.69  kcal  mol"1  with 
HOD  in  which  the  OH  bond  has  been  excited  by  4  quanta.  The  necessary  energy  to  create 
products  is  obtained  mostly  from  excitation  of  the  HOD  reagent  so  that  there  is  sufficient 
available  energy  (24.3  kcal  mol1)  to  excite  up  to  OD(v=3)  and  OH(v=2)  products.  Product 
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internal  energy  distributions  were  probed  with  laser  induced  fluorescence.  Table  IV  shows  the 
experimental  and  theoretical  relative  populations  of  OD(v)  and  OH(v).  The  theoretical  results 
only  use  electronic  state  1.  We  note  that  classical  trajectory  products  can  be  assigned  to  quantum 
bins  that  can  exceed  the  total  energy  available  to  products,  a  problem  inherent  in  conventional 
histogram  binning  in  QCT.  As  described  in  Ref.  20,  most  of  the  initial  vibrational  energy  in  the 
OH  bond  of  HOD  is  retained  in  the  OH  product,  and  the  OD  functions  like  a  spectator  in  the 
reaction.  The  theoretical  results  agree  well  with  the  measured  OD  populations,  but  there  is  only 
qualitative  agreement  with  measurements  for  the  OH(v)  populations.  The  theoretical  results  are 
largest  for  OH(v=2)  but  there  is  appreciable  population  in  v=0  and  v=l. 


TABLE  IV.  Relative  populations  of  OD(v)  and  OH(v)  for  the  measurements  of  Ref.  20. 


OD(Ref.  20) 

OD  (theory) 

OH(Ref.  20) 

OH(theory) 

v=0 

0.93±0.03 

0.81 

0.17 

V=1 

0.07±0.03 

0.16 

0.00±0.04 

0.34 

v=2 

— 

0.03 

1.0±0.04 

0.49 

Figure  14  shows  the  relative  populations  of  OD(v=0,j)  and  OH(v=2,j)  states,  for  the 
measurements  of  Ref.  20  and  the  present  theoretical  results.  There  is  fairly  good  agreement  for 
OD(v=0,j)  populations,  but  there  is  considerable  disagreement  for  the  OH(v=2,j)  populations. 
Due  to  the  histogram  binning  method  in  the  presently  implemented  QCT,  the  theoretical  results 
predict  OH(j>10)  states  which  are  not  energetically  allowed  and  this  is  extremely  problematic  in 
the  present  case  where  the  vibrational  bins  are  quite  wide.  This  problem  can  be  addressed  with 
the  use  of  gaussian  binning,  something  we  plan  to  carry  out  in  future  work.  In  this  case,  Gaussian 
binning  could  have  a  dramatic  impact  on  the  rotational  distributions.  We  also  show  in  Fig.  14  the 
relative  populations  for  OD(v=0,j)  and  OH(v=2,j)  for  a  rotational  temperature  of  700  K  which  is 
fairly  close  to  the  temperature  implied  by  the  relative  collision  energy.  The  700  K  temperature 
populations  are  a  fairly  good  fit  to  the  data,  and  this  could  suggest  that  collisional  relaxation 
plays  a  major  role  in  determining  these  populations. 
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IV.  FUTURE  WORK 


The  present  calculated  results  for  reaction  (1)  are  well  suited  for  planning  and  analyses  of 
hyperthermal  O-atom  measurements  which  should  be  possible  in  the  near  future.  If  final  state 
resolution  is  possible,  such  measurements  would  be  an  ideal  way  to  assess  the  accuracy  of  the 
present  surfaces  and  dynamics  results  and  explore  the  rich  dynamics  of  this  system.  The  present 
surfaces  can  also  be  used  to  investigate  direct  collisional  excitation  of  H20(X  Ai)  by  0(  P), 
reaction  (2).  For  further  analysis  of  the  laboratory  and  space-based  spectral  data,  direct 
collisional  excitation  results  will  be  critical.  This  is  probably  best  carried  out  with  quantum 
scattering  approaches,  as  resolution  of  final  quantum  states  for  polyatomic  systems,  although 

i  -j 

done  previously,  presents  special  problems  with  classical  approaches.  It  is  also  possible  to 

2  2 

model  OH(XTI)  +  0H(X“1T)  collisions,  the  reverse  of  reaction  (1)  with  the  present  surfaces, 

although  the  present  surfaces  only  involve  three  of  the  four  electronic  states  on  the  OH  +  OH 

2  2 

side.  Since  OH(XTI)  +  OH(X  n)  collisions  to  make  H2O  are  exothermic,  measurements  of  the 
reverse  of  (1)  are  more  experimentally  accessible.  We  should  mention  also  that  the  present  work 
does  not  consider  singlet  surfaces  nor  electronic  coupling  which  could  be  important  in 
comparing  to  detailed  measurements.  Finally,  we  note  that  the  present  surfaces  form  the  lower 
part  of  a  large  electronic  state  manifold  leading  to  electronically  excited  OH(A,  X  ).  The 
present  approaches  and  surfaces  together  with  information  from  Ref.  12  could  be  used  as  the 
basis  to  construct  global  surfaces  with  couplings  leading  to  electronically  excited  OH  observed 
in  several  measurements.10’11 
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FIGURE  CAPTIONS 

Figure  1.  Energy  level  diagram  of  low-lying  triplets  (solid  lines)  and  singlets  (dashed  lines) 
showing  reagents,  transition  states  and  products  for  O  +  FLO  collisions.  The  transition  state 
geometries  and  energy  levels  are  for  triplets  based  on  the  work  of  Ref.  17.  The  energy  and 
corresponding  relative  collision  velocities  are  shown  in  parenthesrs  with  respect  to  the  0(  P)  + 
FLOIX,1  A|)  asymptote  at  zero  energy.  The  energies  of  the  fundamental  vibrational  modes  of 
FLO  are  shown  for  reference.  In  the  present  study,  we  follow  the  three  lowest  triplet  states  of 
adiabatically  0(3P)  +  H20(X,1Ai)  to  OH(2fI)  +  OH(2fI). 


3  I 

Figure  2.  Internal  coordinates  used  in  the  description  of  the  0(  P)  +  H20(X,  Ai)  potential  energy 
surfaces.  Active  coordinates  (dashed  lines)  are  the  set  {x,y,z}  and  spectator  coordinates  (solid 
lines)  are  the  set  {u,9,x}.  (a)  illustrates  the  {x,y,z,u,9}  ‘in-plane’  coordinates  and  (b)  illustrates 
the  ‘out-of-plane’  x  coordinate  looking  along  the  Oa-Hb  bond. 
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•3 

Figure  3.  Fitted  potential  energy  surfaces  of  the  three  lowest  triplet  states  of  the  0(  P)  + 
FFOfX,1  A|)  system.  On  each  frame,  the  abscissa  is  the  x-coordinate  and  the  ordinate  is  the  y- 
coordinate.  Each  column  shows  a  different  electronic  state.  Each  row  shows  a  different  fixed  z- 
coordinate  (80°,  140°  and  240°).  The  spectator  coordinates  {u,9,x}  are  fixed  at  their  transition 
state  values  as  given  in  Table  I.  Contour  levels  are  at  {1,5,10,15,20,25,30,40,50}  kcal  mol'1. 
Distances  are  in  a.u. 


Figure  4.  The  t  dependence  of  the  active  potential,  Va  for  the  three  lowest  triplet  states  of  the 

3  1 

0(  P)  +  H20(X,  Ai)  system.  Flere  it  is  defined  with  respect  to  the  energy  of  Va  for  state  1  at  the 
transition  state  given  in  Table  I. 


Figure  5.  The  V$  portion  of  the  spectator  potential,  Vs,  as  a  function  of  0  for  the  three  lowest 
triplet  states  of  the  0(3P)  +  FFCKX/Ai)  system: - state  1; - ,  state  2;  — ,  state  3. 

Figure  6.  Full  potential,  Vsym,  for  planar  geometries  of  the  0(3P)  +  FFCXX,1  A})  system  where  the 
reagent  H2O  bond  lengths  and  angles  are  fixed  at  the  transition  state  geometry,  as  shown  in  Table 
I  for  the  three  lowest  triplet  adiabatic  states.  The  position  of  the  reagent  O-atom  is  varied  along 
the  abscissa  and  ordinate  coordinates.  Distances  are  in  a.u.  Contour  values  are  at 
{1,5,10,15,20,30,40,  50}  +  16.8  kcal  mol1  (which  is  the  value  of  the  potential  for  the  O-atom  far 
away  for  state  1). 

3  1 

Figure  7.  Rate  constants  as  a  function  of  temperature  for  the  reaction,  0(  P)  +  H20(X,1A1)  -> 
OH(X,2n)  +  OH(X,2n).  |_|-  experimental  results  of  Ref.  22;  X,  experimental  results  of  Ref.  23; 
O,  total  calculated  results  including  contributions  from  all  three  states;  O,  contribution  from 
state  1;  V,  contribution  from  state  2;  A,  contribution  from  state  3. 


Figure  8.  Total  calculated  OFI  cross  section  as  a  function  of  collision  energy  for  the  reaction 
0(3P)  +  H20(X,1Ai)  — »  OH(X,2n)  +  OH(X,2n)  according  to  Eq.  12:  O,  total  calculated  results 
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including  contributions  from  all  three  states;  <0,  contribution  from  state  1;  V,  contribution  from 
state  2;  A,  contribution  from  state  3.  Relative  collision  velocity  is  shown  on  the  top  scale. 


Figure  9.  (a)  Calculated  vibrationally  resolved  cross  sections  as  a  function  of  collision  energy  for 
the  reaction  O  +  H20  — >  OFI(v)  +  OFI(v)  according  to  Eq.  12.  For  a  given  v,  contributions  from 
both  OFI(v)  products  are  summed,  (a),  O,  results  at  5  km  s'1  collision  velocity;  A,  results  at  7 
km  s'1  collision  velocity;  A,  results  at  8  km  s'1  collision  velocity;  V,  results  at  9  km  s'1  collision 
velocity;  □,  results  at  1 1  km  s'1  collision  velocity,  (b),  same  as  (a)  except  results  only  include 
OH(v)  products  formed  from  the  incoming  O-atom  reagent  (active  OFI).  (c),  same  as  (a)  except 
results  only  include  OH(v)  products  formed  from  reagent  H2O  O-atom  (spectator  OH). 


Figure  10.  Calculated  rotationally  resolved  cross  sections  as  a  function  of  collision  energy  for  the 
reaction  O  +  H20  — »  OH(j)  +  OH(j)  according  to  Eq.  12.  For  a  given  OH  rotational  level  j, 
contributions  from  both  active  and  spectator  OH  products  are  summed  over  all  OH(v)  vibrational 
levels.  O,  results  at  5  km  s'1  collision  velocity;  A,  results  at  7  km  s'1  collision  velocity;  A, 
results  at  8  km  s'1  collision  velocity;  V,  results  at  9  km  s'1  collision  velocity;  [_|,  results  at  1 1 
km  s'1  collision  velocity. 


Figure  1 1  .Calculated  relative  populations  at  8  km  s'1  relative  collision  velocity  for  the  reaction  O 
+  H20  -»  OHactiveG)  +  OHspcctatorO).  — □ — ■,  presents  results  for  OHactiveG);  -  -  -  surprisal 

theory,  Eq.  16,  with  <9, =0.0;  — O — ,  presents  results  for  OHspectatorG);  -  -  O-  -  surprisal  theory, 
Eq.  16,  with  0r=  -0.41. 


Figure  12.  Calculated  product  energy  fraction  as  a  function  of  collision  energy  for  the  reaction 
0(3P)  +  H20(X,'Ai)  — >  OH(X,2n)  +  OH(X,2n).  (a)  total  summed  over  OH  active  and  spectator 
contributions,  — A — ,  translational  energy  fraction,  fj;  — O — ,  rotational  energy  fraction,  fr; 
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vibrational  energy  fraction,  fv.  (b)  active  OH  and  spectator  OH  contributions  to  the 
internal  energy  fraction. 


Figure  13.  Calculated  differential  angular  cross  sections  as  a  function  of  collision  velocity.  O, 
results  at  5  km  s'1  collision  velocity;  A,  results  at  7  km  s'1  collision  velocity;  A,  results  at  8  km 
s'1  collision  velocity;  V,  results  at  9  km  s'1  collision  velocity;  □,  results  at  1 1  km  s'1  collision 
velocity.  To  be  consistent  with  previous  results,  we  follow  the  convention  for  Eq.  12  for  the  total 
cross  section. 


Figure  14.  Rotational  relative  populations  for  the  reaction  0(3P)  +  HOD(4vOH)  — >  OD(v,j)  + 
OH(v,j).  (a)  calculated  results  for  OD(v=0)  product,  -  -  Q-  -,  experimental  results  of  Ref.  20; 
— O — ,  present  calculated  results;  —A—,  population  of  OD(v=0,j)  at  700  K.  (b)  same  as  (a) 
except  for  OH(v=2,j). 
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Fig.  1 
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Fig.  3 
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Fig.  4 
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Fig.  5 
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Fig.  6 
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Fig.  7 
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Fig.  8 
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Fig.  10 


41 

'Approved  for  public  release;  distribution  unlimited. 


Relative  Population 


11 


42 

'Approved  for  public  release;  distribution  unlimited. 


Energy  Fraction 


Fig.  12 
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Fig.  13 
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Fig.  14 
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